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I. INTRODUCTION 

Surveys of type la supernovae and detailed mappings of t 
provide strong evidence that the universe is accelerating 



re cosmic microwave background 
l|. To explain this acceleration 



within the framework of Einstein's theory of General Relativity requires the existence of 

n 

some form of "dark energy" which has positive energy density and negative pressure |2|]. 

One common model for dark energy Q is to treat it as a pervasive, homogenous per- 
fect fluid with equation of state p = wp. Cosmic acceleration demands that w < —1/3, 
and observations from the Wilkinson Microwave Anisotropy Probe in conjunction with su- 
pernova surveys and baryon acoustic oscillations measurements place the current value at 
w = — 0.999^056 0| . Although this is consistent with the effective equation of state for a 
cosmological constant, w = —1, we cannot rule out the possibility that our Universe contains 
"phantom energy" , for which w < — 1 . 

If w is a constant and w < — 1 then general relativity predicts that as the universe expands 
the phantom energy density increases with the result that in a finite amount of proper time 
the phantom energy density will become infinite and the universe will expand by an infinite 
amount. All bound objects, from clusters of galaxies to atomic nuclei, will become unbound 
as the Universe approaches this future singularity, aptly called the "Big Rip" (J). 

A simple model of a spacetime with a Big Rip singularity can be obtained by considering 
a spatially flat Robert son- Walker spacetime, with metric 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ). (1.1) 

At late times the phantom energy density is much larger than the energy density of classical 
matter and radiation and the solution to the classical Einstein equations is 

a(t) « ai(t r -t)- a , (1.2) 

with 

a = _2/(3 + 3w) > . (1.3) 

Here a\ is a constant and t r is the time that the Big Rip singularity occurs. The phantom 
energy density is 

P P h = ^-a 2 (t r -ty 2 . (1.4) 
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Note that both the scale factor and the phantom energy density become infinite at t — t r . 

In addition to the singularity in the model described above, other classes of models con- 
taining somewhat milder phantomlike future singularities have been identified. Barrow |5[ 
constructed a class of models called "sudden singularities" in which w is allowed to vary as 
w oc (t r — i) a ~ l for < a < 1. In these models, the pressure and scalar curvature diverge 
at time t r but the energy density and scale factor remain finite. Other types of singular be- 
havior were found in Refs. js _ S|- Nojiri, Odintsov and Tsujikawa [9] came up with a general 
classification scheme for future singularities. The strongest singularities are classified as type 

I and Big Rip singularities fall into this class. The sudden singularities are examples of type 

II singularities. Two other classes, type III and type IV were also identified. For type III 
singularities both the energy density p and the pressure p diverge at the time t r but the scale 
factor a remains finite. For type IV singularities, the scale factor remains finite, the energy 
density p and the pressure p go to zero at the time t r , but divergences in higher derivatives of 
H = a/a occur. Other classification schemes for cosmological singularities have been given 



in Refs. 



10|-ll2|. 



At times close to t r in each of the above cases it is possible that quantum effects could 
become large and that the backreaction of such effects could moderate or remove the final 
singularity. One way to investigate whether this would occur is to compute the energy density 
for the quantum fields in the background geometry of a spacetime with a final singularity. 
Then a comparison can be made between the phantom energy density and the energy density 
of the quantized fields. A second way is to solve the semiclassical backreaction equations to 
directly see what effects the quantum fields have. 

on 

Nojiri and Odintsov p, H3[ studied the backreaction of conformally invariant scalar fields 
in the cases of sudden singularities and Big Rip singularities and found that quantum effects 
could delay, weaken or possibly remove the singularity at late times. In Ref. |9J Nojiri, 
Odintsov and Tsujikawa used a model for the dark energy with an adjustable equation of 
state to find examples of spacetimes with type I, II, and III singularities. They then solved 
the semiclassical backreaction equations and found that the singularities were usually either 
moderated or removed by quantum effects. 

Calderon and Hiscock 14j investigated the effects of conformally invariant scalar, spinor, 
and vector fields on Big Rip singularities by computing the stress-energy of the quantized 
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fields in spacetimes with constant values of w. Their results depend on the value of w and 
on the values of the renormalization parameters for the fields. For values of w that are 
realistic for our universe they found that quantum effects serve to strengthen the singularity. 
Calderon 15[ made a similar computation in spacetimes with sudden singularities and found 



that whether the singularity is strengthened or weakened depends on the sign of one of the 
renormalization parameters. 

Barrow, Batista, Fabris and Houndjo [lij considered models with sudden singularities 
when a massless, minimally coupled scalar field is present. They found that in the limit 
t — > t r , the energy density of this field remains small in comparison with the phantom 
energy density. Thus quantum effects are never important in this case. Batista, Fabris 

n 

and Houndjo [17| investigated the effects of particle production when a massless minimally 
coupled scalar field is present in spacetimes where if is a constant. To do so they used a state 
for which Bunch and Davies 18| had previously computed the stress-energy tensor. They 
found that the energy density of the created particles never dominates over the phantom 
energy density. 

Pavlov 19( computed both the number density of created particles and the stress-energy 
tensor for a conformally coupled massive scalar field for the case in which w = — 5/3. It was 
found that backreaction effects are not important for masses much smaller than the Planck 
mass and times which are early enough that the time until the Big Rip occurs is greater 
than the Planck time. 

In this paper we compute the energy densities of both massless and massive scalar fields 
with conformal and minimal couplings to the scalar curvature in spacetimes with Big Rip 
singularities in which the parameter w is a constant. While our calculations are for scalar 
fields, it is worth noting that both massive and massless conformally coupled scalar fields can 
be used to model spin 1/2 and spin 1 fields, and in homogeneous and isotropic spacetimes 



the massless minimally coupled scalar field can serve as a model for gravitons [20l . 121 ]. 

For conformally invariant fields the natural choice of vacuum state in homogeneous and 
isotropic spacetimes is the conformal vacuum 22| ■ For all other fields there is usually no nat- 
ural choice. However, it is possible to define a class of states called adiabatic vacuum states 
which, when the universe is expanding slowly, can serve as reasonable vacuum states 221 ] . 
They can be obtained using a WKB approximation for the mode functions, and they are 
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specified by the order of the WKB approximation. It has been shown that the renormalized 
stress-energy tensor for a quantum field is always finite if a fourth order or higher adiabatic 
vacuum state is chosen. 1 

Here we compute the renormalized energy densities of the quantum fields in fourth order 
or higher adiabatic states and investigate their behavior as the universe expands. One focus 
is on the differences that occur for the same field in different states. We find in all cases 
considered that the asymptotic behavior of the energy density is always the same for fields 
with the same coupling to the scalar curvature, regardless of whether they are massless or 
massive and regardless of what states the fields are in. Fields with minimal coupling to the 
scalar curvature have a different asymptotic behavior than those with conformal coupling. 

We also address the question of whether and under what conditions the energy density of 
the quantized fields becomes comparable to the phantom energy density. We find that for 
fields in realistic states for which the energy density of the quantized fields is small compared 
to that of the phantom energy density at early times, and for spacetimes with realistic values 
of w, there is no evidence that quantum effects become large enough to significantly affect 
the expansion of the spacetime until the spacetime curvature is of the order of the Planck 
scale or larger, at which point the semiclassical approximation breaks down. 

In Sec. II the quantization of a scalar field in a spatially flat Robertson- Walker spacetime 
is reviewed along with a method of constructing adiabatic states. In Sec. Ill, the energy 
density for massless scalar fields with conformal and minimal coupling to the scalar curvature 
is discussed and a comparison is made with the phantom energy density. For the massless 
minimally coupled scalar field a proof is given that one particular state serves as an attractor 
state in the sense that for a large class of states, the energy density of the field asymptotically 
approaches the energy density it would have if it was in the attractor state. In Sec. IV 
numerical calculations of the energy density for massive scalar fields with conformal and 
minimal coupling to the scalar curvature are discussed. A comparison is made with both the 
phantom energy density and the energy density of the corresponding massless scalar field. 
Our main results are summarized and discussed in Sec. V. Throughout units are used such 

1 In this paper we generalize the definition of an n'th order adiabatic vacuum state to include all states 
whose high momentum modes are specified by an n'th order WKB approximation but whose other modes 
can be specified in any way. 
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that h—c — G—1 and our sign conventions are those of Misner, Thorne, and Wheeler 23]. 



II. QUANTIZATION IN SPATIALLY FLAT ROBERTSON WALKER 

SPACETIMES 

We consider a scalar field obeying the wave equation 

(□ - m 2 - £R)(j) = 0, 



(2.1) 



where m is the mass of the field and £ is its coupling to the scalar curvature R. It is 
convenient for our calculations to use the conformal time variable 

r* di 



V 



(2.2) 



We expand the field in terms of modes in the usual way 22] 

1 f d 3 k 



a(v) J (2ir)\ 

with the creation and annihilation operators satisfying the commutation relations 

[a k , a k /] = [a k , a k ,] = 
[a k ,al] = a(k-k') . 

The time dependent part of the mode function satisfies the equation 



(2.3) 



6 / a 



with a constraint given by the Wronskian condition 



(2.4) 



(2.5) 



(2.6) 



Throughout primes denote derivatives with respect to the conformal time 77. 

We restrict our attention to states for which the stress-energy tensor is homogeneous 
and isotropic. Thus the stress-energy tensor is uniquely specified by the energy density 
(p) f , = — (T °) and the trace (T). The unrenormalized energy density is y 



(P)u 



47r 2 a 4 J 



dkk'U |Vfc| 2 + 



1\ a 



12 



\ipk 



(2.7) 
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Renormalization is to be accomplished through the use of adiabatic regularization 25k28|. 
Following the prescription given in Ref. 29| the renormalized energy density is 



(P) r = (P)u - (P)d + (P)c 



(2.8a) 



with 2 



and 



(P)c 



47r 2 a 4 



oo r | 

dkk 2 <k + 



4ir 2 a ., x 



dkk 2 { — 
k 6 



m 2 a 2 



4 4 

ma 



1\ 3a 



121 



6 / a 2 



1\ 3m 2 a' 2 
{ -6 — 



+ If - i) 



2880vr 2 

4 

m 
~64vr 2 

,2 



288tt 



/~r o 



1 , //xV 
2 +1 ° g W 



(1) #o° 



6 7 288^ 



+ 



m 



167r ; 



/-> 



(1) #n° 







3 + log 



32tt 2 



\i 2 a 2 
4A 2 " 



2 + log 



/x 2 a 2 
4A 2 " 



9 a! 2 a" 
4tt 2 a 7 



36a'"a' 72a"a' 2 18a" 2 

5— + = + =— 



3 a 



/ 1 



3a 



12 



(2.8b) 



(2.8c) 



(2.8d) 
(2.8e) 
(2.8f) 



Note that the value of (p) is independent of the arbitrary cutoff A. For a massive field fi = m 
whereas for a massless field \x is an arbitrary constant. The renormalization procedure for 
the trace (T) is given in Ref. 29[ and follows similar lines. For the purposes of this paper, 
we are concerned only with the energy density. 



2 Note that Eq.(9a) of Ref. 29( has a misprint. The term on the third line which is proportional to m 2 
should be multiplied by a factor of (£ — 1/6). 
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A. WKB Approximation and Adiabatic States 



A formal solution to the time dependent part of the mode equation can be given using a 
WKB expansion. First make the change of variables 



i>k{r)) 



: exp 



'/ 



-/ / W(r])dr] 

"no 



where r] is some arbitrary constant. Substituting into Eq. (12 .5ft gives 

1\ a" 1 fW" 3W' 2 \ 



(2.9) 



W 



k' z + m 2 a 2 + 6 U - 



(2.10) 



6 J a 2\W 2W 2 J ' 
Approximate solutions to the above equation may be obtained using an iterative scheme, 
where each successive iteration is given by 



W, 



(<<) 



(n-2) 



3 
2 



W, 



(n-2) 



(2.11) 



with 



W, 



(o) 



(k 2 + m 2 a 2 ) 2 . 



(2.12) 



Taken in conjunction with Eq. (12. 9p . this can be used to define an nth-order WKB approx- 
imation to the mode equation, 



1 



exp 



I Wi n \ri)dri 
'vo 



(2.13) 



<2Wt\v) 

In general, this approximation is valid when k is large and/or the scale factor varies slowly 
with respect to rj. In these cases it will differ from the true solution by terms of higher order 
than n. 

Exact solutions of the mode equation can be specified by fixing the values of ^fc and ip' k at 
some initial time 7] in such a way that the Wronskian condition (12.61) is satisfied. One way 
to do this is to use an nth-order WKB approximation to generate these values. For example 
if (I2.13P is used then 



(2.14) 
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Note that, while it is acceptable to include derivatives of in the expression for tjj' k , it is 
not necessary because such terms are of order n + 1. 

If Eq. (I2.14p is substituted into Eq. (12.71) and then the quantity (To ) is computed, it 
is found that there are state dependent ultraviolet divergences in general for zeroth-order 
and second-order adiabatic states. There are no such state dependent ultraviolet diver- 
gences for fourth- or higher-order adiabatic states. This is true for all times because in 
a Robertson- Walker spacetime an nth-order adiabatic state always remains an nth-order 



adiabatic state 



221. 



In what follows we are interested in investigating a large class of states. Thus as mentioned 
in footnote 1 in the Introduction we generalize the definition of an nth-order adiabatic 
vacuum state to include all states whose high momentum modes are specified in the above 
way by an nth-order WKB approximation but whose other modes can be specified in any way. 
However, we shall not consider states in this paper which result in an infrared divergence of 
the energy density. 

III. MASSLESS SCALAR FIELDS 

A. The Massless Conformally-Coupled Scalar Field 

Setting £ = 1/6 and m = in Eq. (12.51) gives the mode equation for the conformally 
coupled scalar field 

< + A:Vfc = (3.1) 

which has solutions 

p —ikr\ pikrj 



The conformal vacuum is specified by the state with a = 1 and /3 = 0. Following the 
renormalization procedure given in Eq. (I2.8al) . we note that the contribution from the modes 
(p) u is exactly canceled by the counter-terms in (p) d . Thus, the renormalized energy density 
is given by 

<">'= <»- = 28^ (-i<W + »«/'). (3.3) 
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To determine the point at which quantum effects become important one can compare the 
energy density of the quantum field with the phantom energy density. This can be done 
analytically by using Eqs. (I3.9P and (I3.10P to write Eq. (I3.3P in terms of the scalar curvature 
as 

/ \ 1 27w 2 + 18w - 5 d2 

{p) r = 34560^ (1 - 3wy R ■ (3 ' 4) 
The phantom energy density in terms of the scalar curvature is given by 

Pph = 8n(l-3w) R ■ (3 ' 5) 
Since R diverges at the time t r it is clear that eventually the energy density of the conformally 
invariant field becomes comparable to and then larger than the phantom energy density. The 
two densities are equal at the time when 

1 — 3w 

R = 4320tt . 3.6 

27w 2 + 18w - 5 V ; 

For the semiclassical approximation to be valid we need the spacetime curvature to be much 

less than the Planck scale or R -C 1. From Eq. fl3.6p one finds the two densities are equal 

when R — 1 if 

w = - - (3 + 2160vr + 2V6 + 64 8077 + 1166400tt 2 ) w -1500 . (3.7) 

Thus for any physically realistic value of w the energy density of the scalar field will remain 
small compared to the phantom energy density until the scalar curvature is above the Planck 
scale. For example, for w = —1.25, the energy density of the scalar field equals the phantom 
energy density at R ~ 4400. 

B. The Massless Minimally-Coupled Scalar Field 

Setting m = £ = in (12. 5p gives the mode equation for the massless minimally coupled 
scalar field, 

€ + (k 2 " 7) V* = 0. (3.8) 
In terms of the conformal time rj the scale factor is 

a( V ) = a Q {-rj)~\ (3.9) 
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with 



7 = ~ 



l + 3w 



> 



and ao a positive constant. 

The general solution to this equation is given in terms of Hankel functions, 



MV) = *-P- (aH^ il/2) (-k V ) + PH^ (1/2) (-k V ) 
Substitution into the Wronskian condition (12. 6 j) gives the constraint 



(3.10) 



(3-11) 



lal 2 - 



(3.12) 



A special case is the solution with a = 1 and (3 = 0. The stress-energy tensor was 
computed analytically for the state specified by this solution by Bunch and Davies jl8 ]. 
Therefore, in what follows, we will refer to this state as the Bunch-Davies state. Bunch and 
Davies found that 3 



(P) 



BD 



1 



2880vr 2 

1 



(1) #o° + (3) #o° 



+ 



11527T 2 

1 



log ( 4 ) + V>(2 + 7) + i>0- ~ 7) + t 



[-24DR + 2ARR ° + 3R 2 ] 



R 



13824tt 2 

1 351w 2 + 54w - 65 



967r 2 a 2 ?7 2 



+ 



69120vr 2 

1 + D2 



256tt 2 (3w - 1) 



(3w - l) 2 
R 



R 2 



Qw 



4' 

+ 3 



.(3.13) 



i + 3^y vi + sw, 

The mass scale \x in the above expression is an arbitrary constant and ip is the digamma 
function. 

As in the case for the massless conformally coupled scalar field, one can compare Eq. ( 13. 13f) 
with Eq. (13. 5p and see that eventually the energy density of the scalar field will become 
comparable to and then greater than the phantom energy density. A graphical analysis 



3 In the process of our calculation, we discovered a misprint in Eq. (3.36) of Ref. |lSj] which is repeated in 
Eq. (7.52) of Ref. 22j. The final term in both equations should be multiplied by a factor of 2. 
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shows that for w > —140 the two energy densities become comparable at a time when R > 1. 
For example, if w = —1.25 and fi = 1, the energy density of the scalar field equals the 
phantom energy density at R w 180. 

The Bunch Davies state turns out to be an attractor state in the sense that the energy 
density for all homogeneous and isotropic fourth order or higher adiabatic states approaches 
(p)f D . This type of behavior was found in Ref. [30] in de Sitter space for scalar fields in 
the Bunch-Davies state and we use the same type of argument here as was used in that 
paper to establish it. First note that since the renormalization counter terms are the same 
for any choice of a and (3 in Eq. (13. lip , the renormalized energy density can be written in 
terms of the energy density of the Bunch-Davies state plus remainder terms. By substituting 
Eq. (13. lip into Eq. (12.71) . setting m = £ = 0, and using the constraint given in Eq. (I3.12p . 
one finds that for this class of states the energy density can be written as 



(P)r = (P)r° + m 



(3.14) 



with 



■sir a* 



dkk 4 



r(l) 



+ Re (a/3* (V\ {-kr,)f + (tf^(-^)) 2 ) } . (3.15) 



The integral in the above expression can be split into three parts, 



Ifa) = h{v) + h{v) + hiv) 



= t4 f X dkk 4 {...} (3.16a) 



Sira 4 



h(v) = ir4 / dkk 4 {...} (3.16b) 



57ra 4 j x 



-77) 



sira 4 



h(v) = / dkk 4 {...} (3.16c) 

l-Z/v 



where A and Z are positive constants. To have a fourth order adiabatic state, |/3| must fall 
off faster than ci/k A for large enough values of k and any positive constant c±, thus A is 
chosen such that \(3\ < 1/k 4 when k > A. The constant Z is chosen such that < Z « 1. 
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To evaluate the first integral, we use the series expansion 



x 



71 



(l)' v + o(x-^ 2 ) + o(xn 



(3.17) 



This is valid when x = —krj < Z . We note that at late enough times Z > —Xi] for any choice 
of Z and A. Using the expansion (13.171) in Eq. ( 13.151) . one finds to leading order that 



h(v) 



^ [r( 7 + |)] : 



-v) 



2-1 



dk {k 3 ~ 2 "< [\/3f 



-Re (a/3*)} + O [(-k V ) 2 } + O [(-kr,) w ] } . 



(3.18) 



For states with no infrared divergences the integral is finite, so we find that Ii{ff) — > as 
(-77) -> 0. 

For the second integral, due to the time dependence in the upper limit cutoff, we find 
that all orders of the expansion given in Eq. (I3.17P contribute, so a different treatment is 
necessary. Using the fact that for k > A we have |/3| 2 < 1/k 8 and \at/3*\ < c 2 /k 4 for some 
positive constant c 2 , we can put the following upper bound on the integral in Eq. (I3.16bj) . 



\Hv)\ < 



■v) 



-z/v 



lira 4 j x 



+ (n+? 



dk 



-krj) 



k 4 



+ c 2 



J, 



-krf) 



7+5 
2~ 



-krj) 



-krf) 



+2c 2 



J^(-kr,)Y^(-kr]) + J^(-k V )Y^(-kr,) 



7+ 



(3.19) 



Assuming 7 7^ |, one can use the standard series expansions for Bessel functions to show 
that the terms in the integrand are all of the form 



-77) 



47+1 



-z/v 



dk ( ^ + c 2 



n,n'=0 



\n+n' 



-krf) 



(a 1 +a 2 )p+2n+2n' 



n\n'\ T(aip + n + l)T(a 2 p + n' + 1) 



(3.20) 



where «i and a 2 independently take on the values ±1, p takes on the values 7 ± n> an d ci 
and c 2 are constants. After integration, Eq. (I3.20p becomes 
-\) n+n ' 1 



n,n'=0 



n!n'! r(«ip + n + l)T(a 2 p + n' + 1) 



x 



-rf) ci 



+ 



c 2 



Z 4 ((«! + a 2 )p + 2n + 2ri - 3) (cki + a 2 )p + In + 2n' + 1 
c~i 



2-(ai+a 2 )p+2n+2n'+l 



A 4 ((ai + a 2 )p + 2n + In' - 3) 



c 2 



(«i + a 2 )p + 2n + 2n' + 1 



(ai+Q 2 )p+2n+2n.'+l 



(3.21) 
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Note that the above series converges. The smallest power of (—77) is (— 7/) 27 and occurs when 
a i — 0C2 = — 1 and p = 7 + |. It comes from the n = n' = term. Thus, Ii{rj) — >■ as 

(-v) -> 0. 

In the case 7 = ~, we must use the appropriate series expansions for Yq(— krj) and Yi(—krj), 
so Eq. (I3.2ip will contain different terms. However the argument is similar and the result 
that ^(t/) — )■ as (—77) — > is the same. 

For the third integral, using \(3\ 2 < l/k s and \oe/3*\ < C2/A; 4 , and changing variables to 
x — —krj gives the bound 

|l3(r7)| < 



-v) 



47+4 



z x4 



H {1 \(x) 

7-2 



7+5' 



C 2 (-7/) 



47 



(i) 



X . 



+ fi- 



rm 



x . 



(3.22) 



By noting that 



Hv (x) ~ 1/ ^fx at large x, it is clear that these integrals are finite. Further, 



each integral is independent of (—77), so the contribution from this piece goes like some 
constant times (— r/) 47 . Thus, l{r\) — > as {—rf) — > and this completes our proof that 
{p)r° ^ s an attractor state. 

To demonstrate this behavior numerically, we compare the energy density from Eq. (13. 13[) 
which occurs for a = 1 and /3 = to the renormalized energy density of a different fourth 
order adiabatic state. One way to specify a fourth order adiabatic state is to use a fourth 
order WKB approximation as discussed in Sec. Ill Al However, for the spacetimes we are 
considering, using a fourth order WKB approximation in Eqs. (I2.13P and f |2.14j) in the limit 
k — > results in an infrared divergence of (p). To avoid such a divergence, we use a zeroth 
order WKB approximation for the modes with < k < ko , for some positive constant ko. 
According to the generalized definition we give in Sec. I, these are still fourth order adiabatic 
states. A given state is specified by the value of ko and the time t]q at which the matching 
in Eq. fl2.14l) is done. 

As shown in Fig. [TJ the energy density of the state specified by choosing k = 0.2 and 
7/0 = —100, approaches that of the Bunch Davies state as expected. A number of other 
numerical calculations were done for various values of fco and 7/0 with the same qualitative 
result. Note that for all the numerical calculations that were done for this paper, w = —1.25 
and a « 22.6. 
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FIG. 1: The dashed line represents the phantom energy density p p h in the case that w = —1.25. The 
two solid lines represent the energy density (p) r for two states of the massless minimally coupled 
scalar field. The central line is the Bunch Davies attractor state. The oscillating line is a fourth 
order adiabatic state with ko = 0.2 and 770 = —100. For both states the mass scale ll = 1 has been 
used. 



It is of course possible to construct states for which, at times close to r] , the energy 
density is large compared to the phantom energy density. However, there is no reason to 
expect that the energy density of quantized fields at early times should be comparable to 
that of the phantom energy density since that is not the case today. Thus our numerical 
results provide evidence that if the initial energy density is small compared to the phantom 
energy density at times close to 770 then, for realistic values of w, it will remain small during 
the period when the semiclassical approximation should be valid. 
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IV. MASSIVE SCALAR FIELDS 



A. Massive Conformally Coupled Scalar Fields 



Setting £ = 1/6 in Eq. (12. 5p gives the mode equation for a massive conformally coupled 
scalar field, 



In this case, we do not have an analytic solution for the mode equation and there is no 
obvious choice for the vacuum state. Instead, we choose our initial state using the method 
outlined near the end of Section IIIIBI and compute solutions to the mode equation and the 
energy density numerically. What we find is that there is an initial contribution to the energy 
density from the modes of the field in addition to the contribution from Eq. (13. 3p ; however, 
as the Universe expands this contribution redshifts away and the energy density approaches 
the energy density of the massless conformally coupled scalar field. An example is shown in 
Fig. [2] for a field with mass m = 0.005 in a state specified by ko = 0.1 and 770 = —100. We 
have tested this for several different initial states and different masses and found consistent 
behavior in all cases. 

As with the massless minimally coupled scalar field, it is possible to construct states with 
a large initial contribution to the energy density and these states can be ruled out by the 
same argument as given in Section IIIIBI Our numerical results provide evidence that if the 
energy density of the quantized field is small compared to the phantom energy density at 
times close to rjQ, then, for realistic values of w, it will remain small during the period when 
the semiclassical approximation should be valid. 



The last case we consider is that of massive, minimally coupled scalar fields, which obey 
the mode equation 



Again, we choose our initial state using the method outlined near the end of Section IIIIBI 
and evolve the modes forward in time numerically. 



< + (k 2 + m 2 a 2 ) i) k = . 




B. Massive Minimally-Coupled Scalar Fields 




(4.2) 
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FIG. 2: The dashed line represents the phantom energy density p p h in the case that w = —1.25. 
The upper solid line is (p) for a massive {m = 0.005) conformally coupled scalar field in the 
fourth order adiabatic state with k$ = 0.1 and r/o = —100; the lower one is {p) r for the massless 
conformally coupled scalar field. 



Qualitatively, the behavior here is similar to the behavior found for the massless, mini- 
mally coupled field in a fourth order adiabatic vacuum state. At late times, the contribution 
from the mass terms is small compared to other terms and the energy density approaches 
that of the massless minimally coupled scalar field in the Bunch-Davies vacuum state. We 
note that there is an arbitrary mass scale p, present in the renormalized energy density (13.131) 
for the massless scalar field; the convergence shown in Fig. [3] occurs when p, is set equal to 
the mass of the massive field. At very late times the value of p, does not affect the leading 
order behavior of (p) for the massless minimally coupled scalar field and thus, for any value 
of pi, the energy density of the massive field asymptotically approaches that of the massless 
field in all cases considered. 
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FIG. 3: The dashed line represents the phantom energy density p p ^ in the case that w = —1.25. 
The central solid line is the energy density (p) r for a massless minimally coupled scalar field in 
the Bunch Davies state with mass scale p, = 0.001. The oscillating line is the energy density for a 
massive (m = 0.001) minimally coupled scalar field in the fourth order adiabatic state with ko = 1 
and 770 = — 50. Note that the quantity \p\ rather than p has been plotted since the energy density of 
the massive field is negative at 77 = 770, oscillates between positive and negative values for a period 
of time, and then becomes positive definite at late times. Note that, except at the initial time, the 
times where the curve is nearly vertical when it reaches a lower limit are times at which the energy 
density goes through zero, which is — 00 on the scale of this plot. Thus the curve should really 
extend down to 00 at these points. 

V. SUMMARY AND CONCLUSIONS 

We have computed the energy densities of both massless and massive quantized scalar 
fields with conformal and minimal coupling to the scalar curvature in spacetimes with big 
rip singularities. We have restricted attention to states which result in a stress-energy tensor 
which is homogeneous and isotropic and free of ultraviolet and infrared divergences. For the 
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numerical computations we have further restricted attention to states for which, near the 
initial time of the calculation, the energy density of the quantum field is much less than that 
of the phantom field. 

For the massless minimally coupled scalar field we have shown that the energy density for 
the field in any fourth order or higher adiabatic state for which the stress-energy tensor is 
homogeneous, isotropic, and free of infrared divergences, always asymptotically approaches 
the energy density which this field has in the Bunch-Davies state. In this sense the Bunch- 
Davies state is an attractor state. 

For massive minimally coupled scalar fields numerical computations have been made of 
the energy density for different fourth order adiabatic states and in every case considered 
the energy density approaches that of the Bunch-Davies state for the massless minimally 
coupled scalar field at late times. For conformally coupled massive scalar fields numerical 
computations have also been made of the energy density for different fourth order adiabatic 
states. In each case considered the energy density asymptotically approaches that of the 
massless conformally coupled scalar field in the conformal vacuum state. Thus it appears 
that the asymptotic behavior of the energy density of a quantized scalar field in a spacetime 
with a big rip singularity depends only upon the coupling of the field to the scalar curvature 
and not upon the mass of the field or which state it is in, at least within the class of states 
we are considering. 

Analytic expressions for the energy densities of both the massless conformally coupled 
scalar field in the conformal vacuum state and the massless minimally coupled scalar field 
in the Bunch-Davies state, in spacetimes with big rip singularities have been previously 



obtained 
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22] and are shown in Eqs. ( 13. 3 j) and (13. 13D . To investigate the question of 



whether backreaction effects are important in these cases, the energy density of the scalar 
field can be compared to the phantom energy density to see if there is any time at which 
they are equal or at least comparable. Then one can determine whether the semiclassical 
approximation is likely to be valid at this time by evaluating the scalar curvature and seeing 
whether or not it is well below the Planck scale. We have done this and find that for the 
conformally coupled field the two energy densities are equal at the point when the scalar 
curvature is at the Planck scale if w ~ —1500. For the minimally coupled scalar field this 
combination occurs if w ~ —140. Thus for w <C —1500 for the conformally coupled field 
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and w <C —140 for the minimally coupled field, one expects that backreaction effects may 
be important at times when the scalar curvature is well below the Planck scale. However, 
the values of w which satisfy these constraints are completely ruled out by cosmological 
observations. 

Another way in which the energy density of a quantum field can be comparable to the 
phantom energy density at scales well below the Planck scale is to construct a state for 
which this is true. There is no doubt that such states exist. The analytic and numerical 
evidence we have is that over long periods of time the energy density of a conformally or 
minimally coupled scalar field in such a state would decrease and at late enough times 
become comparable to that of a massless field in the conformal or Bunch-Davies vacuum 
state respectively. More importantly one can ask whether any such states exist which are 
realistic for the universe that we live in. This seems unlikely because today is certainly 
an early time compared to t T if our universe does have a big rip singularity in its future; 
however the energy densities of quantized fields today are much less than that of the dark 
energy. Thus it would be necessary for the energy density of a quantum field to be small 
compared to the phantom energy density today and then to grow fast enough to become 
comparable to it well before the Planck scale is reached. This type of behavior seems highly 
artificial, particularly since it is not what happens for a massless scalar field in the conformal 
or Bunch-Davies vacuum states. Thus we find no evidence which would lead us to believe 
that backreaction effects due to quantum fields would remove a big rip singularity in our 
universe if, indeed, such a singularity lies in our future. 
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